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Abstract 

We discuss the interplay of screening, absorption and regeneration effects, on the quantum mechanical evolution of 
quarkonia states, within a time-dependent harmonic oscillator (THO) model with complex oscillator strength. We 
compare the results with data for /?aa/^aa(CNM) from CERN and RHIC experiments. In the absence of a measure- 
ment of cold nuclear matter (CNM) effects at LHC we estimate their role and interpret the recent data from the ALICE 
experiment. We also discuss the temperature dependence of the real and imaginary parts of the oscillator frequency 
which stand for screening and absorption/regeneration, respectively. We point out that a structure in the J/i/' sup- 
pression pattern for In-In collisions at SPS is possibly related to the recently found X(3872) state in the charmonium 
spectrum. Theoretical support for this hypothesis comes from the cluster expansion of the plasma Hamiltonian for 
heavy quarkonia in a strongly correlated medium. 
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1. Introduction 

The J/i/r suppression at CERN and RHIC has been 
quantified by measurements of the nuclear modification 
factor Raa \M- This set of experimental data is com- 
monly normalized to the baseline of cold nuclear matter 
(CNM) effects /?aa(CNM). The results are scaled with 
the multiphcity of charged particles ^j^^g llllllll. 

A key problem to be solved in heavy-ion collisions 
is that the usage of heavy quarkonia states as a probe 
for the diagnostics of the quark-gluon plasma (QGP) re- 
quires the knowledge of the baseline for their produc- 
tion and evolution characteristics in situations when a 
QGP is absent. For a review see, e.g., |j6j] and references 
therein. Current issues in the experiment-theory dia- 
logue are summarized, e.g., in [7;]. Aspects ofquarko- 
nium production at LHC are discussed, e.g., in islsi and 
new information for its 7?aa(CNM) baseline are given in 

The suppression of charmonium (more general, 
heavy quarkonium suppression), due to color screen- 
ing, was suggested to be a signal of QGP formation 



1 1 1LI12I1 . The works lll3Lll4lll5ll have studied the disso 



ciation of quarkonium (charmonium and bottomonium) 



by screening in a hot plasma. The binding energy of 
the heavy quark-antiquark pair in the quarkonium state 
vanishes at the respective Mott temperature where the 
bound state merges the continuum of scattering states 
Olg, 11711 . Below the Mott temperatures, the binding 
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energies are already sufficiently lowered so that colli- 
sions with particles from the medium may have suffi- 
cient thermal energy to overcome the threshold for im- 
pact dissociation of quarkonium iHEH. Both ef- 
fects tend to wash out the pattern of sequential sup- 
pression for heavy quarkonia states expected from the 
"classical" picture of the Mott effect 12111 . The descrip- 
tion of quarkonia states in a hot QGP medium in the 
vicinity of the critical temperature should therefore treat 
bound and scattering states on an equal footing. This 
is appropriately achieved within a thermodynamical T- 
matrix approach, which has been developed to address 
the spectral properties of quarkonia ll22ll as well as open 
flavor meson states 11231 124 12511 . It is also worth not- 
ing that inelastic collisions are responsible for absorp- 
tion of J/i/' in a hadron gas of light n and p mesons since 
the corresponding cross sections can be sufficienly large 
I 20l, 26 , 27 , 28 , 29ll . Once the partial densities of open 
charm mesons becomes nonnegligible, also the reverse 
processes of charmonium regeneration in channels like, 
e.g., D + D ^ + p will occur. J/i// regeneration shall 
give important contributions to charmonium production 
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in heavy-ion collisions already at RHIC and the more at 

the LHc iIMlIi Ills SIM Sill 



Another very interesting alternative is the channel 
D + D* ^ J/iff + p II14II where as a "molecular" bound 
state of D and D* with binding energy below 1 MeV the 



exotic state X(3872) is discussed ||7 ll37ll38ll which was 
recently discovered by BaBar ll39ll and subsequently 
confirmed by other experiments. 

It has been shown that as long as medium effects 
can be embodied in a gaussian action, the QQ propa- 
gator obeys a closed temporal evolution equation whose 
large-time behavior is governed by an effective poten- 
tial. The latter, besides screenin g, displ ays also an imag- 
inary part related to collisions ll40[ kill . The determina- 
tion of this imaginary part also provides insight into how 
the anisotr opy of the medium relates with the width of 



OEV 



the states 11421 14311 . The use of complex potentials for 



describing data of charmonium suppression in pp and 
pA collisions has been documented in 1 44 , 45 , 4a, 47 1 . 
An application of a time-dependent harmonic oscillator 
(THO) model with a complex oscillator frequency to ex- 
perimental data of charmonium production in AA colli- 
sions has recently been given |3 48]. This application 
predicted the formation of resonances nearby the onset 
of charmonium suppression (Mott transition) ll49ll . This 
model exhibits confinement in the sense that the poten- 
tial grows quadratically with the separation distance of 
the QQ pair. 

In this work we consider such quantum-mechanical 
model for studying the role of static screening and ab- 
sorption/regeneration kinetics for the quarkonium time 
evolution. The THO model for the production of QQ 
pairs is a generalization of the Matsui formula [i50il to a 
complex QQ potential. The quantum mechanical evo- 
lution, as a function of the time t is encrypted in the 
evolution amplitude U{r,ro,t) for the initial QQ state 
ipQQ. We have examined two alternatives for the initial 
state: it can be taken as a delta function or as 

a Gaussian ||49|1 . The time dependence originates from 
the temperature evolution of the medium T{t) which is 
assumed to obey total entropy conservation. In order 
to account for hot medium effects produced in AA colli- 
sions, the temperature dependence of the complex oscil- 
lator frequency aP'iT) stands for screening, absorption 
and regeneration in a unified approach. By taking differ- 
ent ansatze for u^{T), we show that the THO model can 
describe the behavior of the measured J/i/r yield over the 
expected CNM effects i?AA/^AA(CNM), as a function of 
the multiplicity of charged particles. Our method can be 
applied for more realistic cases, even for a QQ Hamil- 
tonian containing a complex confined potential in better 
agreement with Lattice QCD (LQCD) as the given in 



|42l|43|]. 

We pay special attention to the In-In data 101 taken 
at CERN SPS since they have the highest precision. 
Interestingly, they exhibit a peculiar shape ("wiggle") 
lITii I49I1 . The wiggle turns out to be important since 
it appears around the estimated temperature for disso- 
ciation of charmonium (~ 1.3 with Tc - 0.154 
GeV f^JS), which can be defined as the temperature for 
its Mott transition. Our novel result with the complex 
THO model suggests that such peculiar behavior can be 
mapped into 3[(y^(r)], in other words into the imagi- 
nary part of the QQ potential. We examine whether the 
production of the X(3872) resonance can contribute to 
the in-medium QQ hamiltonian. 

2. Suppression factor for quarkonium 

The survival probability (suppression factor) for l/tfr 
is defined in i50ll 



(1) 



This expression was published in bra-ket representation 
by Cugnon et al. [45]. This is a generaUzation of the 
Matsui approach to an arbitrary initial state ll49ll . The 
intermediate QQ state is given by 

^QQir.t)^ j flfVo0e2(ro,O)f/(^,^o,O, (2) 

where <pQQ{fo, 0) represents the initial QQ state at the 
position ro and time f = 0. In the following we shall 
assume that the quantum number n, on the quarkonium 
state i//„, stands for the principal, angular and azimuthal 
quantum numbers. The suppression factor can be gen- 
eralized to the case of nonsingular complex potentials 
and be applied to the QGP diagnostics in collisions of 
heavy ions with mass number A when identified with 
the experimentally determined quantity 



7?aa(CNM) 



(3) 



where 7?aa(CNM) accounts for the CNM effects 
from charmonium absorption in cold nuclear mat- 
ter and modification of charm production by shadow- 
ing/antishadowing of gluon distribution functions in the 
CM of the colliding nuclei. Both effects are accessi- 
ble by analysis of pA collision experiments, see IH 0] 
and references therein. We restrict our discussion here 
to ground state charmonium at rest in the QGP medium 
iPT = y = 0) so that the discussion of Lorentz boost 
effects on the formation process can be omitted. We 



2 



include, however, the feed-down from higher charmo- 
nia states which are discussed more in detail elsewhere 

HEi. 

The next section is dedicated to the application of the 
generalized Matsui approach. To this end, we start cal- 
culating the evolution amplitude U{r,ro,t). As an ex- 
ample we consider the THO Hamiltonian to model the 
QQ interaction. 

3. Time-dependent harmonic oscillator model 



methods ||48[|49[|58| 



3/2 



, (7) 



2msm{(f>{t)- (piO)) 

where the classical action is exactly calculated by 
= I ■ r - ro ■ roj lUgllsl]. Additioanlly the initial 



Sci = - r- n) ■ foj 

contitions for the Ermakov equation ([6} are defined as 



P(0) = 



MP) 



,P(0) = 0. 



For our discussion of the quantum mechanical evolu- 
tion of quarkonia in an evolving QCD plasma state, we 
will employ here a generalization of the harmonic oscil- 
lator model ifsoll to a time-dependent one with nonsingu- 
lar complex squared oscillator frequency (THO model). 
The merit of such a model is its simplicity and trans- 



4911 . Aspects 



parency as well as its tractability 11141 148 , 
of an optical potential for the propagation of charmonia 
through a medium have been already discussed, e.g., for 



cold nuclear matter in ||46, 52 1 and for a quark-gluon 
plasma in 1451 15311 . We consider the time-dependent 
Hamiltonian for heavy quarkonia in the form 



H{t) = ^ + ^w^(T)r2(T) , 



(4) 



where ji - mQ/2 is the reduced mass and nig the heavy 
quark mass. The squared complex oscillator frequency 
aj^(T) has an implicit time dependence due to the 
temperature evolution T{t) of the system surrounding 
the evolving heavy quarkonium state. The quadratic 
dependence of the imaginary part of the (optical) 
oscillator potential is motivated by the phenomenon of 
color transparency, see also 115 411 and references therein. 



The general classical trajectories for the Hamiltonian 
© can be found in Q El [H IH S 1131 and refer- 
ences therein. It is given by a linear combinations of the 
two solutions 



ri.iir) = p(t) exp(±i0(T)) , 0(t) 



(5) 



The amplitude p(t) fulfills the Ermakov equation 14 

5^ El 



p(T) + C0\t) p{T) - 



= 0, 



(6) 



for which exact solutions exist. The initial and final con- 
ditions of the oscillator are defined as r(0) - ro and 
r{t) - r, respectively. This allows to evaluate the transi- 
tion amplitude as a function of time using path integral 



4. Applications to Heavy Ion Collisions 

In the following we show that the anomalous J/i^ sup- 
pression in SPS, RHIC, and LHC experiments can be 
simultaneously described with the natural assumption 
that above the critical temperature the relevant screen- 
ing of the quarkonium interaction can be parameter- 
ized with a temperature dependent, complex oscillator 
strength. The time evolution of the temperature itself 
will be given by longitudinally boost invariant (Bjorken 
scaling) hydrodynamic evolution of a fireball volume 
Vije) under entropy conservation 



s(T{TeW{T,) = const ; V{t,) = Arz(T,) 



(8) 



The temperature dependence of the entropy density s{T) 
is taken from recent lattice QCD simulations p ill wh ich 
are well parameterized by the simple ansatz lll4[ |49 1 



s{T) = 9.0 
see Fig.[T] 



1 + tanh 
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Lattice data: Borsanyi et al., arxiv:1 007.2580 
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Figure 1: Recent results (data points) for the equation of state from 
lattice QCD jslll compai'ed to the fit fomiula ((9) employed here. 
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Figure 2: Initial temperature Tq and time interval t = Tf — tq for the 
evolution of the QQ state from fomiation at tq to freeze-out at Tf as 
a function of the measured charged particle yields per pseudorapitity 
interval (dNch/dr]),po. 



The initial values of temperature correspond to initial 
entropy densities given by the Bjorken formula (|59 1 



^0 = 



3.6 



dA^\ 



(10) 



where the transverse overlap area At is taken in this 
work as |fl4i|49ll 



At = (4.8 fmO 



1=0 



(11) 



We found the parametrization by fitting the average val- 
ues of the overlap area for PbPb. The latter was obtained 
from Monte Carlo calculations performed in Ref. Q. 
Additionally the transformation iHlslI - l-04x 

I X ( ^) has been used. 

^ V d'; /rj=0 



We are interested in the time evolution of the quarko- 
nium state in the cooling fireball medium described by 
the thermal history ^(t) starting from t - tq - 1.0 
fm/c until the freeze-out of the effective T HO p otential 
parameters at Tf = T{Tf) = 0.193 GeV iH. Here, 
we defined t/ = TQs{T())/s{Tf). In Fig.|2]we show the 
initial temperature Tq - T{to) and the duration of the 
quarkonium evolution until its freeze-out t - Ty -tq as a 
function of the charged particle multiplicity per pseudo- 
rapidity interval idNch/dri)r,=o characterizing the initial 
conditions of the fireball evolutiorQ- 



' In using Eq. (T) for describing anomalous suppression due to 
a plasma medium switched on at a time tq after the creation of the 
QQ pair implies that the evolution of the latter before tq is governed 
by a hermitian Hamiltonian which does not yet suffer a modification. 



We consider the most suitable parametrization of the 
THO model with complex squared frequency given bjH 
o/iT) = %[oj^(T)] + i 3[ci/(T)]. In this case screen- 
ing, absorption and regeneration can be described si- 
multaneously. Screening occurs when the real part of 
the squared frequency is ^R[w-(f)] < w^, where o)^ is 
the oscillator frequency for the J/if/ in vacuum. The 
case of absorption appears when the im aginary part is 



3[6/(0]< 0. The regeneration scenario |l30l |3Il |32l |33| 



corresponds to 3[to^(T)] > in the THO 
model which stems from the coupling of the charmo- 
nium hamiltonian to open charm channels. In addition 
we shall consider a gaussian wave function for initial 
QQ state. Details are presented in ll49ll . 

The ground state of quarkonium can be identified 



with the IS state of the harmonic oscillator ||48 



with - 



given by lAooo(r) = i/'ooo(0)exp^j^j 
Calculations of the suppression factor with Dirac delta- 
shaped initial QQ state were performed in lll4ll on the 
basis of Ref. 14811 . In the following we use results 
for the J/i/r survival probability obtained in ll49ll us- 
ing a gaussian shaped initial state and discuss the re- 
lation between parametrizations of the temperature de- 
pendent complex squared oscillator strength to experi- 
mental data on anomalous J/^A suppression from SPS, 
RHIC and LHC. In particular, we discuss how account- 
ing or not for the "wiggle" in the In-In data at SPS when 
parameterizing uP'iT) will change predictions of the sur- 
vival probability in the LHC domain. 



5. Anomalous suppression from SPS to LHC 

Anomalous suppression, the deviation from unity of 
experimental data for the J/i/r production ratio normal- 
ized to the expectation accounting for CNM effects (O, 
is considered as a key indicator for QGP formation 
in heavy-ion collisions. This effect, first observed at 
CERN SPS for Pb-Pb collisions at yfs ^ 17 GeV, has 
been qualitatively confirmed by RHIC experiments with 
Au-Au interactions at Vi = 200 GeV whereby three 
particularly interesting observations were made 

(i) the suppression is stronger at forward and back- 
ward rapidities rather than at midrapidity where the 
particle densities are highest, 

(ii) the onset of anomalous suppression and its depen- 
dence on centrality scales with the charged particle 



^Note that in Ref. 
quency has been used. 



a different convention for the THO fre- 
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density at midrapidity rather than with energy den- 
sity. 

(iii) the rather precise data of the NA60 collaboration 
for In-In collisions show a dip in the centrality de- 
pendence of the anomalous J/i/r suppression ratio 

la. 

While (i) is caused mainly by antishadowing and to 
some extent by geometry II60II . the second finding is 
still not understood. For the puzzling In-In dip a sug- 
gestion has been made in ll4ll where this feature could 
be reproduced within the generalized Matsui approach 
by a subtle interplay of screening and absorption in the 
parametrization of the temperature dependence of the 
oscillator frequency. For the details see ll49l . Accord- 
ing to this picture, the dip reflects a nonmonotonous 
temperature behavior of the confining potential due to 
a resonance-like contribution in %[a)^{T)]. This calcu- 
lation is rather susceptible to changes in the absorptive 
part in the complex oscillator strength, basically as a 
consequence of the Dirac delta distribution for the ini- 
tial QQ state. 

In the present work we investigate the role of 
3 [oP'iT)] in reproducing the dip in In-In data, see Fig.|3] 
and Fig. |4] The parametrization used here allows a de- 
scription of screening, absorption and regeneration si- 
multaneously. Additionally, we employ for the initial 
QQ state the gaussian shape given in |49]. For the LHC 
data on J/ifr suppression we employ recent results of the 
ALICE collaboration I16IJ. For estimating the CNM ef- 
fects which are not measured yet at LHC, we deduce an 
error band guided by theoretical estimates in Ref. 1162 1 
adn shown as the cyan hatched region in Figs.[3]|5] 

For an additional discussion of the behavior of the 
imaginary part of the QQ potential shown in Fig. |4]we 
refer to a separate analysis of the QQ plasma Hamilto- 
nian l63ll and the discussion of the role of the X(3872) 
particle and the spectral function of p meson in a quark 
meson plasma in this context Il4ti49i1 . It is worth noting 
that the wiggle pattern in the In-In data is located in the 
region of temperatures 1.3 < T/T^ < 1.4 Ollii]. The 
influence of this wiggle on the charmonium suppression 
in the LHC region with charged particle densities ex- 
ceeding the ones reached at RHIC is explored in Fig.|5] 
In this figure the two LHC predictions are completely 
different above dNchldr] ~ 700 although the complex 
oscillator strength is the same above 0.23 GeV. How 
can this wiggle affect the LHC predictions so strongly? 
One explanation can be seen in the survival probabil- 
ity ([TJ which encodes the quantum time evolution of the 
quarkonium state. As we have described above, the time 
evolution of the complex oscillator frequency is implied 
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Figure 3: Predictions of the quantum meclianical THO model for char- 
monium suppression at CERN and RHIC. The LHC data are described 
by a purely imaginary part of the QQ Hamiltonian. 



by its assumed temperature dependence using the tem- 
perature evolution of the plasma T{t) which follows 
from entropy conservation under Bjorken expansion (|8) 
with the fit of the lattice QCD data for the entropy den- 
sity (in. 

The predictions for LHC are affected by the value 
of the oscillator strength at T - Tj ~ 1.3 where 
the complex frequency oj(To) enters the transition am- 
plitude U{r,ro,t) (see Eq. (|7])) through the condition 
p(0) = 1 / Vw(0) and is determined by the equality 
ojiTo) = a)(Tf). 

6. Conclusion 

Recent studies of heavy quarkonia correlators and 
spectral functions at finite temperature in lattice QCD 
and systematic T-matrix approaches using QCD moti- 
vated finite-temperature potentials support that heavy 
quarkonia dissociation shall occur in the temperature 
range L2 < Td/T^ < 1.5 whereby the interplay of 
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Figure 4: The same as Fig.|3] but fitting the "wiggle" in the In-In data 
by varying the imaginary part of the frequency. The dip is mostly de- 
scribed by absorption and regeneration on the region from 0. 193 GeV 
to 0.21 GeV. The screening is rather visible above Tf = 1.3 T^, as a 
consequence the production of charmonium decreases despite of the 
presence of the imaginary part. In order to exhibit a regeneration pat- 
tern in the domain of LHC beyond the RHIC data, the imaginary pait 
of the squared oscillator frequency for temperatures above 0.23 GeV 
shall be larger than the corresponding one given in Fig.0 



both screening and absorption processes is important. 
We have discussed these effects on the quantum me- 
chanical evolution of quarkonia states within a time- 
dependent harmonic oscillator model with complex os- 
cillator strength and compared the results for the sur- 
vival probability with data for /?aa/^aa(CNM) from 
SPS, RHIC and LHC experiments. Besides the tradi- 
tional interpretation, with a threshold for the onset of 
anomalous suppression by screening and dissociation 
kinetics at dNch/drj ~ 300, we suggest an alternative 
arising from the attempt to model the dip of the suppres- 
sion pattern of the rather precise NA60 data from In-In 
collisions at dNch/drj ~ 150 - 250. We suggest that this 
dip indicates the true threshold for the onset of anoma- 
lous suppression due to the coupling of charmonium to 
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Figure 5: Influence of the wiggle on the charmonium suppression at 
the LHC region. This region correspond to temperatures above 0.23 
GeV see reference |49|| . The two LHC predictions are completely dif- 
ferent above dNch/drj ~ 700 although the complex oscillator strength 
above 0.23 GeV is the same. The difference is due to the different 
values of 3[w(Ty )]. 



the channel with the recently discovered X(3872) 

state. Although some details have been worked out in 



038114911 . the theoretical basis for supporting this hypoth- 
esis has apparently been developed in plasma physics 
with the concept of a plasma Hamiltonian for nonrel- 
ativistic bound states like heavy quarkonia when they 
are immersed in a medium dominated by strong corre- 
lations like bound states. In the context of the theory 
of strongly correlated plasmas the formation and disso- 
ciation of bound states can be systematically addressed 
within cluster expansion techniques. With the advent of 
the new quality of the data from the ALICE experiment 
at LHC which provide yields of hidden {QQ = J/i/r, T, 
...) and open heavy flavor (Qq = D, B, ...) mesons, it 
may be possible for the first time to measure the heavy 
quark association degree nQg/nggn^Q. Such quantity 
is of advantage for the QGP diagnostics not only be- 
cause it is known for plasma diagnostics in stellar atmo- 
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spheres, semiconductor plasmas and nuclear matter, but 
in particular since it provides a new baseline for quanti- 
fying the J/i/r suppression effect. We will return to this 
issue in a separate work. 

In the present work we have illustrated an alterna- 
tive to the models of charmonium kinetics and statisti- 
cal hadronization in the description of data for anoma- 
lous J/i^ suppression that utilizes the quantum mechani- 
cal formulation of quarkonia production with a complex 
Hamiltonian. To illustrate the method, a simple THO 
model with complex oscillator frequencies has been pre- 
sented whereby the temperature dependence of the real 
and imaginary parts resulting from fitting heavy-ion col- 
lisions data of SPS, RHIC and LHC experiments on 
quarkonium production may allow to discuss the inter- 
play of screening, absorption and regeneration of heavy 
quarkonia in these experiments. It is obvious that this 
tool may be further refined by using realistic tempera- 
ture dependent complex potentials as motivated by lat- 
tice QCD studies. The quantum mechanical approach 
can also be exploited, e.g., to study the relative abun- 
dances of different excited states (Saha equation aspect 
of plasma diagnostics) and to study provide insights 
about the role of resonances at the open charm thresh- 
olds (most prominently the X(3872)) for the regenera- 
tion of J/i^ in an analogy to the role of the Hoyle state in 
element synthesis. 
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